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1  Objectives 

'  l 

The  main  objective  of  this  work  has  been  the  development  and  maturation  of  adjoint-based, 
viscous  design  techniques  for  Aerodynamic  Shape  Optimization  (ASO)  of  complete  aircraft 
configurations.  We  have  followed  a  systematic  approach  in  the  development,  validation  and 
testing  of  our  methods.  In  addition,  we  have  carried  out  some  preliminary  work  in  the 
development  of  &  mathematical  environment  for  high-fidelity  aero-structural  optimization  as 
a  first  step  towards  the  realization  of  a  high-fidelity  multidisciplinary  optimization  capability. 
A  number  of  sigfiificant  milestones  achieved  in  this  process  are  detailed  in  this  final  report. 
Details  of  the  wbdc  carried  out  can  be  found  in  a  number  of  conference  papers  and  archival 
publications  that  are  quoted  in  the  References  section  [1J. 

\ 

2  Summary 

The  effort  to  de|4lop  a  framework  for  viscous  aerodynamic  shape  optimization  of  complete 
aircraft  configurations  has  yielded  a  number  of  important  contributions  during  the  length 
of  this  grant.  The  following  sections  reflect  these  contributions.  We  start  with  a  thorough 
description  of  the:  adjoint  equations  for  the  compressible  Reynolds-Averaged  Navier-Stokes 
equations  with  gome  mention  of  implementation  details.  Once  the  first  implementation  of 
the  viscous  adjoint  was  completed,  we  carried  out  a  thorough  accuracy  study  of  the  gradients 
produced  by  the!  method  with  very  encouraging  results.  In  addition,  we  also  spent  much  of 
the  second  year  establishing  the  importance  of  continuous  vs.  discrete  adjoint  formulations. 
Our  conclusions-  are  detailed  below  and  establish  the  validity  of  the  continuous  approach, 

j  \ 
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where  we  use  thje:  subscripts  I  and  II  to  distinguish  between  the  contributions  associated 
with  the  variation  of  the  flow  solution  Sw  and  those  associated  with  the  metric  variations 
SS.  Thus  [MJ\j  and  \Pj\t  represent  and  |2  with  the  metrics  fixed,  while  S Mu  and  SVu 
represent  the  contribution  of  the  metric  variations  SS  to  SM  and  6V. 

In  the  steady;  state,  the  constraint  equation  specifies  the  variation  of  the  state  vector  Sw 

by  h 

6R=j-6(Ft-F«)  =  0.  (3) 

Here,  R(w,  S )  is  ijthe  residual  of  the  flow  equations,  and  F(  and  Fvi  are  the  inviscid  and 
viscous  fluxes.  Here,  also,  SR,  SF{  and  SFvi  can  be  split  into  contributions  associated  with 
6w  and  SS  using!  ithe  notation 

i  ;  SR  as  SRj  +  SRjj  (4) 

SFi  =  [Fiw}j6w  +  SFUJ 

,  SFyt  =  [Fviw]j  6w SFyitf.  (5) 

Multiplying  by  a  Co-state  vector  ip  and  integrating  over  the  domain  produces 

JvVSRd^  =  J^iPT^i6(Fi-Fvi)d'Di  =  0.  (6) 

!  *■ 

If  ip  is  differentiable  the  terms  with  subscript  I  may  be  integrated  by  parts  to  give 

<‘  f  iPT8RdV^  +  j  mipTS  {Fi  -  Fvi)j  dBi  (7) 

I!  -  Jv^HF<-F,<hdV'=0.  (8) 

Since  the  left  hanid  expression  equals  zero,  it  may  be  subtracted  from  the  variation  in  the 
cost  function  (1)  to  give 

ti  =  Sin-  jvxpT6Rn-  jB[8M-niipT6{Fi-Fvi)l]dBs 

;;  +  i\s'p+%i  «-■&),]<»«•  (9) 

Now,  since  ip  is  an  arbitrary  differentiable  function,  it  may  be  chosen  in  such  a  way  that  SI 
no  longer  depends  explicitly  on  the  variation  of  the  state  vector  Sw.  The  gradient  of  the 
cost  function  can;  then  be  evaluated  directly  from  the  metric  variations  without  having  to 
recompute  the  variation  Sw  resulting  from  the  perturbation  of  each  design  variable. 

Comparing  equations  (2)  and  (5),  the  variation  Sw  may  be  eliminated  from  (9)  by  equat¬ 
ing  all  field  tenjris  with  subscript  “J”  to  produce  a  differential  adjoint  system  governing 

*  ;;  T 

■ :  [3.  -  a,  J,  +  (P. J,  -  0  to  v.  (10) 

The  corresponding  adjoint  boundary  condition  is  produced  by  equating  the  subscript  “J” 
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boundary  terras  iii  equation  (9)  to  produce 

I’’  "rtF  fa, I-M,],  onB.  (11) 


The  remaining  tprms  from  equation  (9)  then  yield  a  simplified  expression  for  the  variation 
of  the  cost  functjdn  which  defines  the  gradient 

5Iu  +  [  VSRndVt  (12) 

: . 

The  details  tjfi  the  formula  for  the  gradient  depend  on  the  way  in  which  the  boundary 
shape  is  parameterized  as  a  function  of  the  design  variables,  and  the  way  in  which  the  mesh  is 
deformed  as  the  boundary  is  modified.  Using  the  relationship  between  the  mesh  deformation 
and  the  surface  modification,  the  field  integral  is  reduced  to  a  surface  integral  by  integrating 
along  the  coordinate  lines  emanating  from  the  surface.  Thus  the  expression  for  SI  is  finally 
reduced  to 

61=  I  G5SdBs  =  Q6F,  (13) 

j  !.  ^ 

where  7  represents  the  design  variables,  and  Q  is  the  gradient,  which  is  a  function  defined  over 
the  boundary  surface,  The  advantage  is  that  (13)  is  independent  of  5w,  with  the  result  that 
the  gradient  of  /  with  respect  to  an  arbitrary  number  of  design  variables  can  be  determined 
with  only  a  single  flow-field  evaluation  and  a  single  adjoint  evaluation  in  a  design  cycle. 
Then,  the  computational  cost  of  a  single  design  cycle  is  roughly  equivalent  to  the  cost  of  two 
flow  solutions  since  the  adjoint  problem  has  similar  complexity.  When  the  number  of  design 
variables  becomes  large,  the  computational  efficiency  of  the  control  theory  approach  over 
the  traditional  approach  (which  requires  direct  evaluation  of  the  gradients  by  individually 
varying  each  design  variable  and  recomputing  the  flow  field)  becomes  compelling. 

Once  equatibn  (13)  is  established,  an  improvement  can  be  made  with  a  shape  change 

67=-xg, 


where  A  is  positive  and  small  enough  that  the  first  variation  is  an  accurate  estimate  of  SL 
Then  j ; 


si=-xgTg<o. 


After  making  such  a  modification,  the  gradient  can  be  recalculated  and  the  process  repeated 
to  follow  a  path; of  steepest  descent  until  a  minimum  is  reached. 

The  boundary  conditions  satisfied  by  the  flow  equations  restrict  the  form  of  the  left  hand 
side  of  the  adjoint  boundary  condition  (11).  Consequently,  the  boundary  contribution  to  the 
cost  function  M  i.cannot  be  specified  arbitrarily.  Instead,  it  must  be  chosen  from  the  class 
of  functions  which  allow  cancellation  of  all  terms  containing  Sw  in  the  boundary  integral 
of  equation  (9).  j  On  the  other  hand,  there  is  no  such  restriction  on  the  specification  of  the 
field  contribution  to  the  cost  function  V,  since  these  terms  may  always  be  absorbed  into  the 
adjoint  field  equation  (10)  as  source  terms. 

The  detailed;  inviscid  adjoint  terms  have  been  previously  derived  in  [9, 10]  and  a  detailed 
derivation  of  the'  viscous  adjoint  terms  and  the  corresponding  viscous  adjoint  boundary 
conditions  can  ajlso  be  found  in  [2,  3,  11]. 
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4  Accuracy  and  Cost  Study  of  Viscous  Adjoint-Based 
Gradients 

A  continuous  adjjoint  method  for  Aerodynamic  Shape  Optimization  (ASO)  using  the  com¬ 
pressible  Reynolds- Averaged  Navier-Stokes  (RANS)  equations  and  the  Baldwin-Lomax  tur¬ 
bulence  model  was  implemented  and  tested.  The  resulting  implementation  was  used  to 
determine  the  accuracy  in  the  calculation  of  aerodynamic  gradient  information  for  use  in 
ASO  problems.  The  reader  is  referred  to  previous  work  [12]  for  details  of  the  required  math¬ 
ematical  derivations.  The  accuracy  of  the  resulting  derivative  information  is  investigated  by 
direct  comparison;  with  finite-difference  gradients.  The  advantages  of  the  use  of  an  adjoint 
method  over  traditional  finite-differencing  become  apparent  because  of  the  strict  require¬ 
ments  that  the  finite  difference  method  imposes  on  the  level  of  flow  solver  convergence  and 
the  sensitivity  of  the  value  of  the  gradients  with  respect  to  the  choice  of  step  size.  A  paral¬ 
lel  implementation  using  a  domain  decomposition  approach  and  the  MPI  (Message  Passing 
Interface)  standard  is  used  to  reduce  the  computational  cost  of  automatic  design  involving 
viscous  flows.  ;  ’■ 

As  a  result  of  this  study,  the  following  conclusions  were  reached: 

•  The  flow  solver  convergence  requirements  for  the  finite  difference  method  to  produce 
accurate  gradient  information  are  quite  severe  and  substantially  increase  the  compu¬ 
tational  cost  of  this  method. 

I  ' 

•  As  opposed  to  what  may  be  expected  from  direct  truncation  error  analysis,  the  choice 
of  step  size  for  the  calculation  of  derivatives  using  the  finite  difference  method  can 
cause  a  large  error  in  the  approximation  of  the  gradients.  Smaller  step  sizes  (even 
disregarding  the  effect  of  roundoff  error)  are  not  necessarily  better;  the  error  in  the 
resulting  gradients  depends  on  both  the  convergence  level  of  the  flow  solver  and  the 
step  size  cljiosen. 

•  Gradient  information  obtained  using  the  adjoint  method  is  much  less  dependent  on 
the  level  of  convergence  of  the  flow  solver.  Typically  a  convergence  level  2  orders  of 
magnitude  larger  than  in  the  finite  difference  method  can  be  tolerated. 

•  Gradient  information  obtained  using  the  adjoint  method  is  insensitive  to  the  step  size 
chosen  in  the  deformation  of  the  aerodynamic  configuration. 

•  The  adjoiiit  method  requires  only  modest  levels  of  convergence  (1.5  -2.5  orders)  of  the 
adjoint  solver,  thus  reducing  even  further  the  computational  cost  of  this  procedure, 

•  Accurate  gradient  information  can  be  obtained  for  flow  governed  by  the  RANS  equa¬ 
tions  using  the  procedure  outlined  in  previous  work, 

As  a  typical  example  of  the  results  obtained  during  this  study,  Figure  1  shows  the  sensitivity 
of  an  inverse  design  cost  function  to  the  variations  in  50  parameters  that  modify  the  shape 
of  the  airfoil.  As  can  clearly  be  seen,  the  necessary  level  of  convergence  required  for  accurate 
gradients  using  ia  finite-difference  approach  is  quite  a  bit  higher  than  for  the  adjoint  case. 
Not  only  is  the  cost  of  the  adjoint  procedure  independent  of  the  number  of  design  variables 
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i  ■ 

in  the  problem,  tj>iit  also  ,  the  cost  of  the  flow  solution  is  significantly  reduced  due  to  the  less 


la:  Navier-Stokes  Inverse  Design;  Finite  lb:  Navier-Stokes  Inverse  Design:  Adjoint 
Difference  Gradients  for  Varying  Flow  Solver  Gradients  for  Varying  Flow  Solver 
Convergence.  Convergence. 


5  Comparison  of  the  Continuous  and  Discrete  Adjoint 
Approaches  to  Automatic  Aerodynamic  Optimiza¬ 
tion 

The  objective  of  this  study  was  to  compare  the  continuous  and  discrete  adjoint-based  au¬ 
tomatic  aerodyn’ajmic  optimization  approaches.  Our  intention  was  to  study  the  trade-offs 
between  the  corrifilexity  of  the  discretization  of  the  adjoint  equation  for  both  the  continuous 
and  discrete  approaches,  the  accuracy  of  the  resulting  estimates  of  the  gradients,  and  their 
impact  on  the  (jomputational  cost  to  approach  an  optimum  solution.  For  details  on  the 
derivation  and  implementation  of  the  continuous  adjoint  method,  please  refer  to  the  cited 
bibliography.  We  will  briefly  present  the  complete  formulation  of  the  discrete  adjoint  equa¬ 
tion  and  will  then  investigate  their  differences.  The  similarities  between  the  continuous  and 
discrete  boundary  conditions  are  also  explored.  The  results  demonstrate  two-dimensional 
inverse  pressure  I  design  as  well  as  the  accuracy  of  the  sensitivity  derivatives  obtained  from 
continuous  and  discrete  adjoint-based  equations  when  compared  to  finite-difference  gradi¬ 
ents.  More  details  can  be  found  in  reference  [13]. 

The  discrete ’adjoint  equation  is  obtained  by  applying  the  control  theory  directly  to  the 
set  of  discrete  field  equations.  The  resulting  equation  depends  on  the  type  of  scheme  used 
to  solve  the  floW  equations.  The  following  discussion  uses  a  cell  centered  multigrid  scheme 
with  upwind  biased  blended  first  and  third  order  fluxes  as  the  artificial  dissipation  scheme. 
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A  fall  discretization  of  the  equation  would  involve  discretizing  every  term  that  is  a  function 
of  the  state  vector 


rxx  ny 


<j/ =  */«.  +  ££  'PlJ  (R  Mij  +  D  (wh)  > 


(14) 


where  SIC  is  the  discrete  cost  function,  R(w)  is  the  field  equation,  and  D(w)  is  the  artificial 
dissipation  terra,; 

Tterms  multiplied  by  the  variation  5w^  of  the  discrete  flow  variables  are  collected  and 
the  following  is  the  resulting  discrete  adjoint  equation, 


y^Li  — 

&t 


( fry  \M.\T  _  Ax  [iAlr  t*±hl 

Mid)  2 


_  ( Ay  reilT_Aa:  fia1T>)±=u 

la*V U  U 


-i.i  l**J id)  2 

+(W„.  \&]T.  -Afe.., 


-  (AaVj  [fe]y  -■ [S£]y) 

~  [£]y  ”  A^i,y  fe]f,y)  ^ 

+  (A^i^[^~AXr  ■  ^]T 


-  (A%+*  -  Ayk„i  [fi]  y) 


_  ria]J')^3 

vjj  2 


2 


where 


+  (Aav$  [fc]«  -A^-i  [£]*)  ^ 

+H+i,y  -  *<k-id  +  H,y+i  “ 


H+i,y  =  -  iM  “  c<+*,y^+2j 

+3e<+b'  (Vvn,y  -  M  +  c<_|1y^-i,y 


(15) 


(16) 


is  the  discrete  aidjjoint  artificial  dissipation  term  and  V  is  the  cell  area.  The  dissipation 
coefficients  es  and  e4  are  functions  of  the  flow  variables,  but  to  reduce  complexity  they  are 
treated  as  constants. 

In  the  case  of;  an  inverse  design,  6Ie  is  the  discrete  form  of  the  square  of  the  pressure 
integral  around  the  airfoil  surface.  In  contrast  to  the  continuous  adjoint,  where  the  boundary 
condition  appeals1  as  an  update  to  the  costate  variables  in  the  cell  below  the  wall,  the  discrete 
boundary  condition  appears  as  a  source  term  in  the  adjoint  fluxes,  At  cell  i,  2  the  adjoint 
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equation  is  as  follows, 

y&£ii2  — 

9t 

; ;  a  [“^-1,2  Mv  -  V'l-i.a)  -  AT+1,2  (V'i+i.a  -  iM] 
+  2  [— -Bfl  (^<.3  "  ^<,2)]  +  $1 

where  V  is  the  c'ell  area,  $  is  the  source  terra  for  inverse  design, 

i  ij 

;  $  ~  -  (p  -  Pr)A«<)  <fo,2) 

and  s  i1 


(17) 


<i2  =  , 


a/ 


flu; 


i,2 


Aim, 


dg 


dw 


i,  2 


All  the  terms  in  Equation  (17)  except  for  the  source  term  are  scaled  by  the  square  of  A®, 
Therefore,  as  the; mesh  width  is  reduced,  the  terms  within  parenthesis  in  the  source  term 
divided  by  A s<  must  approach  zero  as  the  solution  reaches  a  steady  state.  One  then  recovers 
the  continuous  adjoint  boundary  condition  described  in  earlier  works. 

If  a  first  order*' artificial  dissipation  equation  is  used,  then  equation  (16)  would  reduce  to 
the  term  associated  with  e2.  In  such  a  case,  the  discrete  adjoint  equations  are  completely 
independent  of  the  costate  variables  in  the  cells  below  the  wall.  However,  if  we  use  the 
blended  first  and'  third  order  dissipation,  then  these  values  are  required.  In  practice,  a 
simple  zeroth  ortfer  extrapolation  across  the  wall  produces  good  results. 

Replacing  the  inverse  design  boundary  condition  in  equation  (17)  by  the  discrete  form  of 
the  cost  function  results  in  a  discrete  adjoint  equation  for  drag  minimization. 

As  an  example  of  the  results  in  the  study,  an  inverse  design  test  case  is  briefly  shown 
here.  The  target;  pressure  is  first  obtained  using  the  FL083  flow  solver  for  the  NACA  64A410 
airfoil  at  a  flight  condition  of  M  —  0.74  and  a  lift  coefficient  of  Q  —  0.63  on  a  192  x  32  C-grid. 
At  such  a  condition  the  NACA  64A410  produces  a  strong  shock  on  the  upper  surface  of  the 
airfoil,  thus  making  it  an  ideal  test  case  for  the  adjoint  versus  finite  difference  comparison. 

The  gradientifor  the  continuous  and  discrete  adjoint  is  obtained  by  perturbing  each  point 
on  the  airfoil.  Wb  apply  an  implicit  smoothing  technique  to  the  gradient,  before  it  is  used 
to  obtain  a  direction  of  descent  for  each  point  on  the  surface  of  the  airfoil. 

Figures  (1),  (2),  and  (3)  exhibit  the  values  of  the  gradients  obtained  from  the  adjoint 
methods  and  finite  difference  for  various  grid  sizes.  The  circles  denote  values  that  we  obtain 
by  using  the  finite  difference  method.  The  square  represents  the  discrete  adjoint  gradient. 
The  asterisk  represents  the  continuous  adjoint  gradient.  The  gradient  is  obtained  with 
respect  to  variations  in  Hicks-Henne  sine  “bump”  functions  placed  along  the  upper  and 
lower  surface  of  ithe  airfoil  [14,  12]  •  The  figures  only  illustrate  the  values  obtained  from  the 
upper  surface  starting  from  the  leading  edge  on  the  left  and  ending  at  the  trailing  edge  on 
the  right.  ' " 

Figure  (4)  presents  the  effect  of  the  partial  discretization  of  the  flow  solver  to  obtain 
the  discrete  adjoint  equation.  Here  we  obtain  the  finite  difference  gradients  in  the  figure 


A 


p.ll 


OCT  31  '03  03:31PM  ERA  STANFORD 


Table  1:  Li  norm  of  the  Difference  Between  Adjoint  and  Finite  Difference  Gradient 

without  freezing!  the  dissipative  coefficients.  A  small  discrepancy  exists  in  regions  closer  to 
the  leading  edge:gind  around  the  shock. 

Table  1  contains  values  of  the  Li  norm  of  the  difference  between  the  adjoint  and  finite 
difference  gradients.  The  table  illustrates  three  important  facts:  the  difference  between 
the  continuous  44joint  and  finite  difference  gradient  is  slightly  greater  than  that  between 
the  discrete  adjoint  and  finite  difference  gradient;  the  norm  decreases  as  the  mesh  size  is 
increased;  and  the  difference  between  continuous  and  discrete  adjoint  gradients  decreases  as 
the  mesh  size  is  reduced.  The  second  column  depicts  the  difference  between  the  continuous 
adjoint  and  finite  difference  gradient.  The  third  column  depicts  the  difference  between 
the  discrete  adjoint  and  finite  difference  gradients.  The  last  column  depicts  the  difference 
between  the  discrete  adjoint  and  continuous  adjoint.  As  the  mesh  size  increases  the  norms 
decrease  as  expected.  Since  we  derive  the  discrete  adjoint  by  taking  a  variation  of  the  discrete 
flow  equations,  we  expect  it  to  be  consistent  with  the  finite  difference  gradients  and  thus  to 
be  closer  than  the  continuous  adjoint  to  the  finite  difference  gradient.  This  is  confirmed  by 
numerical  results,"  but  the  difference  is  very  small.  As  the  mesh  size  increases,  the  difference 
between  the  continuous  and  discrete  gradients  should  decrease,  and  this  is  reflected  in  the 
last  column  of  table  1. 

In  conclusion,: 

1.  The  continuous  boundary  condition  appears  as  an  update  to  the  costate  values  below 
the  wall  for.  a  cell-centered  scheme,  and  the  discrete  boundary  condition  appears  as  a 
source  term  in  the  cell  above  the  wall.  As  the  mesh  width  is  reduced,  one  recovers  the 
continuous  adjoint  boundary  condition  from  the  discrete  adjoint  boundary  condition. 

2.  Discrete  adjoint  gradients  have  better  agreement  than  continuous  adjoint  gradients 
with  finite:  difference  gradients  as  expected,  but  the  difference  is  generally  small. 

3.  As  the  m$sh  size  increases,  both  the  continuous  adjoint  gradient  and  the  discrete 
adjoint  gradient  approach  the  finite  difference  gradient. 

4.  The  difference  between  the  continuous  and  discrete  gradient  reduces  as  the  mesh  size 
increases.  , 

I 

i  •• 

5.  The  cost  of  deriving  the  discrete  adjoint  is  greater. 

) 

6.  With  our  search  procedure  as  outlined,  the  overall  convergence  of  the  objective  func¬ 
tion  iB  not;  significantly  affected  when  the  discrete  adjoint  gradient  is  used  instead  of 
the  continuous  adjoint  gradient.  Consequently,  we  find  no  particular  benefit  in  using 
the  discrete  adjoint  method,  which  requires  greater  computational  cost.  However,  we 
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Figure  1:  Adjoint 'Versus  Finite  Dif¬ 
ference  Gradient^?. for  Inverse  Design 
of  Korn  to  NACiA  64A410  at  Fixed 
Ct. 

Coarse  Grid  -  96  x  16,  M  =  0.74, 
Q  =-0.63 


Figure  3:  Adjoint.  Versus  Finite  Dif¬ 
ference  Gradieniisifor  Inverse  Design 
of  Korn  to  NACA  64A410  at  Fixed 

a. 

Fine  Grid  -  25)3;  x  64,  M  -  0.74, 
Cl  0.63 


P.12 

10 


Figure  2:  Adjoint  Versus  Finite  Dif¬ 
ference  Gradients  for  Inverse  Design 
of  Korn  to  NACA  64A410  at  Fixed 
Ci. 

Medium  Grid  -  192  x  32, 

M  =  0.74,  Ct  «  0.63 


Figure  4:  Adjoint  Versus  Finite  Dif¬ 
ference  Gradients  for  Inverse  De¬ 
sign  of  Korn  to  NACA  64A410  at 
Fixed  Ci.  Dissipative  Coefficients 
Not  Frozen 

Medium  Grid  -  192  x  32, 

M  =  0.74,  Ci  =  0.63 
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believe  it  beneficial  to  use  the  discrete  adjoint  equation  as  a  guide  for  the  discretization 
of  the  continuous  adjoint  equation. 

6  Industrial  Applications  to  Three-Dimensional  Wing 
Design  j. 

During  the  period  of  this  grant,  the  three-dimensional  design  method  has  been  used  suc¬ 
cessfully  in  som?  ^significant  industrial  applications  at  Boeing  Long  Beach.  In  addition  to 
other  similar  efforts,  the  method  was  used  to  support  the  studies  of  the  Blended-Wing-Body 
concept  for  long  range  transport  aircraft  which  is  under  study  at  Boeing  with  Air  Force  and 
NASA  support.  [The  method  was  used  to  minimize  the  sum  of  the  computed  drag  coefficients 
at  Mach  0.85  at; three  lift  coefficients  in  the  range  from  0.42  to  0.46.  In  order  to  check  the 
results,  Navier-Stpkes  calculations  were  performed  using  another  CFD  code,  CFL3D,  from 
NASA.  These  calculations  verified  an  improvement  in  the  lift  to  drag  ratio  over  a  range  of 
lift  coefficients  from  0.35  to  0.68,  with  an  improvement  in  the  maximum  lift  to  drag  ratio  of 
4  percent.  At  the, same  time,  the  thickness  of  the  outboard  part  of  the  wing  was  significantly 
increased,  which'  jcould  lead  to  benefits  in  reduced  structural  weight  and  increased  fuel  vol¬ 
ume.  The  result  :of  a  similar  calculation,  checked  with  CFL3D,  can  be  seen  in  the  Figure 
below.  The  comparison  is  between  the  original  design  (dashed  lines)  and  the  optimized  one 
(solid  lines),  respiting  in  more  benign  pressure  distributions  and  an  increase  in  lift-todrag 
ratio  of  4  %.  | 

A  second  application  was  made  to  a  high  wing  transport  design  under  study  at  Long 
Beach.  In  this  case  the  viscous  wing  design  code  was  run  for  a  single  design  point  at  Mach 
0.85.  A  significant  drag  reduction  was  achieved  with  a  calculation  using  60  design  cycles,  with 
only  20  iterations!  of  the  flow  solver  and  20  for  the  adjoint  solver  in  each  design  cycle.  The 
total  CPU  time  using  8  processors  was  57613  seconds.  A  well  converged  analysis  of  the  flow 
requires  400  iterations  that  amount  to  7568  seconds  with  8  SGI  R10000  195  Mhz  processors. 
Thus,  the  optimization  was  achieved  at  a  cost  of  about  7.5  flow  solutions.  Again,  a  benefit 
in  the  lift  to  drig  ratio  was  verified  over  a  Mach  number  range  form  0.845  to  0.865,  with 
even  greater  beriefit  at  the  higher  Mach  number  than  at  the  design  Mach  number  of  0.855. 
These  results  cohfirm  the  utility  of  the  optimization  method  for  real  design  applications.  In 
industrial  use,  the  ability  to  achieve  optimizations  with  computational  costs  equivalent  to  10 
or  even  fewer  fiojw  solutions  is  crucial  to  the  viability  of  the  method.  Research  investigations 
which  lead  to  the  ability  to  construct  design  procedures  with  such  low  computational  cost 
will  be  published'.in  [15]. 

Other  applications  to  supersonic  configurations  have  also  been  presented,  together  with 
details  of  the  algorithms,  implementation,  and  parallelization  [16,  17]. 

|  |! 

7  Preliminary  Application  to  Aero-Structural  Optimiza¬ 
tion 

The  work  contained  in  the  attached  papers  [16,  17, 18, 19]  presents  a  new  framework  for  the 
coupled  optimization  of  aero-structural  systems.  The  framework  permits  the  use  of  high- 
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COMPARISON  OF  CHORDWISE  PRESSURE  DISTRIBUTIONS 
BWBDLR6i/SYN88  FIXED  S.L  AND  QMOD=1.0 


Figure  5:  Sample  Viscous  Design  for  the  Blended-Wing-Body  Airplane. 

fidelity  modeling  :of  both  the  aerodynamics  and  the  structures  and  represents  a  first  step  in 
an  effort  towards  the  development  of  a  high-fidelity  multidisciplinary  optimization  capability. 
The  approach  is  based  on  efficient  analysis  methodologies  for  the  solution  of  the  aerodynamics 
and  structures  ^ubproblems,  an  adjoint  solver  to  obtain  aerodynamic  sensitivities,  and  a 
multiprocessor  parallel  implementation.  We  have  placed  a  geometry  database  representing 
the  outer  mold  Imp  (OML)  of  the  configuration  of  interest  at  the  core  of  our  framework.  Using 
this  geometry  description,  the  information  exchange  between  aerodynamics  and  structures 
is  accomplished  through  an  independent  coupling  of  each  discipline  with  the  OML  database. 
The  framework  'permits  the  later  inclusion  of  other  disciplines,  such  as  heat  transfer  and 
radar  signature^',  with  relative  ease.  Specific  results  from  the  coupling  of  a  finite  volume 
flow  solver  for  the  Euler  and  Reynolds  Averaged  Navier-Stokes  equations  with  two  different 
linear  finite  elertuent  structural  models  were  explored.  Care  was  taken  in  the  treatment  of 
the  coupling  of  the  disciplines  such  that  a  consistent  and  conservative  scheme  is  achieved. 
Direct  comparisons  with  wind-tunnel  data  were  presented  to  demonstrate  the  importance  of 
aeroelastic  solutions.  In  addition,  simplified  design  examples  were  presented  to  illustrate  the 
possible  advantages  of  the  new  aero-structural  design  methodology  in  evaluating  trade-offs 
between  aerodynamic  performance  and  structural  weight  for  complete  aircraft  configurations. 
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This  design  environment  has  been  used  to  perform  RANS  aeroelastic  analysis  of  complete 
configuration  flight  and  wind-tunnel  models  with  an  additional  cost  which  is  less  than  10% 
of  the  cost  of  a  traditional  rigid-geometry  CFD  solution.  These  solutions  can  be  used  to 
determine  a  pridri  whether  significant  aeroelastic  corrections  will  or  will  not  be  needed  for 
the  resulting  wifid  tunnel  data. 

In  addition,  simplified  design  cases  have  been  tested  that  include  the  effect  of  aeroelastic 
deformations  in  the  design  process.  These  cases  have  shown  that  our  design  methodology  is 
able  to  predict  the  correct  trades  between  aerodynamic  performance  and  structural  properties 
present  in  these  types  of  wing  design  problems. 

Finally,  a  structural  stress  penalty  function  was  added  to  the  coefficient  of  drag  of  the 
complete  configuration  to  allow  elimination  of  artificial  thickness  constraints  that  are  typ¬ 
ically  imposed  ifr  aerodynamic  shape  optimization  methods.  This  rudimentary  coupling  of 
aerodynamics  and  structures  in  the  design  not  only  eliminates  the  necessity  to  impose  ar¬ 
tificial  constraint's,  but  also  produces  designs  where  trade-offs  between  aerodynamic  and 
structural  performance  are  considered. 

Further  work -will  focus  on  the  continued  development  of  the  proposed  MDO  frame¬ 
work.  Topics  requiring  significant  research  include  sensitivity  analysis,  optimization  strategy, 
Navier-Stokes  based  design,  use  of  commercially  available  CSM  codes,  multipoint  design,  and 
CAD  integration,  - 

i  |  • 

8  Persoiipel  Supported 

'  Funds  from  this  grant  have  been  used  to  support  the  time  of  Profs.  Antony  Jameson  and 

Juan  J.  Alonso,  pind  Mr.  Siva  Kumaian  (Ph.D.  Candidate  in  the  department  of  Aeronautics 
Astronautics)  | 

9  Honors/ Awards 

•  1 

Session  Chairman.  International  Congress  of  Aeronautical  Sciences  (ICAS)  98  Congress, 
Sidney,  Sept  lj-18,  1998.  AJ 

Invited  Lecture,  i  European  Science  Foundation  Conference  on  Applied  Mathematics  for  In¬ 
dustrial  Flow  Problems,  San  Feliu  de  Guixols,  Spain ,  October  1-3,  1998.  AJ 

Invited  Lecturer:.  Michigan  State  University  Short  Course  on  Computational  Fluid  Dynamics 
Modeling  and  Simulation,  June  16-18,  1999.  AJ 

Lecturer.  Fonda&one  CIME  International  Mathematical  Summer  Center,  Course  on  "Com¬ 
putational  Mathematics  Driven  by  Industrial  Applications ”,  Martina  Franca,  Italy,  June 
81-37,  1999.  AJ : 

Plenary  Lecture ti,  Jth  International  Congress  on  Industrial  and  Applied  Mathematics,  IC1AM 
99,  Edinburgh,  Scotland,  July  6-9,  1999.  AJ 

Invited  Lecture,  i Swedish  Program  in  Scientific  Computing  Annual  Meeting.  Linkoping,  Swe - 
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den.  AJ 

Invited  Lecture.  1  &IT  Chautauqua  Computational  Fluid  and  Solid  Mechanics.  Cambridge, 
MA,  June  2000. :  AJ 

Member  of  NRC\Committee  on  Naval  Hydromechanics.  Washington,  DC.  Report  titled:  “An 
Assessment  of  Naval  Hydromechanics  Science  and  Technology ”  issued.  AJ 

Session  Chairman.  International  Congress  of  Aeronautical  Sciences  (ICAS)  2000  Congress, 
Harrogate,  UK,  August  2000.  AJ  &  JJA 

Session  Chairman.  AIAA  Aerospace  Sciences  Meeting  &  Exhibit,  Reno,  NV,  January  2000. 
JJA 

Invited  Lecture.  [Princeton  University  Mechanical  &  Aerospace  Engineering  Seminar  Series, 
Princeton  University,  December  1999.  JJA 

Invited  Panel  Member.  Computational  Aerosciences  Workshop,  NASA  Ames  Research  Cen¬ 
ter,  April  2000.  J \JA 

Invited  Panel  Afember.  High-Performance  Scalability  Workshop,  Santa  Fe,  New  Mexico. 
May  2000.  JJA  I 

Invited  Panel  Member.  The  History  of  CFD,  AIAA  Aerospace  Sciences  Meeting  &  Exhibit, 
Reno,  NV,  January  2000.  AJ 

Invited  Panel  Member.  The  Future  of  CFD.  AIAA  Fluids  2000  Conference,  Denver,  CO, 
June  2000.  AJ  ,  • 

Invited  Lecture,  j  Computational  Algorithms  for  Aerodynamic  Analysis  and  Design.  Fluent, 
Inc.,  Lebanon,  NH,  August  2000.  AJ 

Invited  Lecture.  ‘  Computational  Algorithms  for  Aerodynamic  Analysis  and  Design.  EADS 
Munich,  Germany,  December  2000.  AJ 

Docteur  Honoris  Causa.  University  of  Paris,  December  2000.  AJ 
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